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Abstract 

Recently the authors used the Kadanoff-Baym non-equilibrium Green's func- 
| tion formalism to derive kinetic equation for the non-condensate atoms, in 



conjunction with a consistent generalization of the Gross-Pitaevskii equation 
for the Bose condensate wavefunction. This work was limited to high tem- 
peratures, where the excited atoms could be described by a Hartree-Fock 
particle-like spectrum. We present the generalization of this recent work to 



low temperatures, where the single-particle spectrum is now described by the 



Bogoliubov-Popov approximation. We derive a kinetic equation for the quasi- 
particle distribution function with collision integrals describing scattering be- 

X" 

tween quasiparticles and the condensate atoms. From the general expression 
for the collision integral for the scattering between quasiparticle excitations, 
we find the quasiparticle distribution function corresponding to local equilib- 
rium. This expression includes a quasiparticle chemical potential that controls 
the non-diffusive equilibrium between condensate atoms and the quasiparti- 
cle excitations. We also derive a generalized Gross-Pitaevskii equation for the 
condensate wavefunction that includes the damping effects due to collisions 
between atoms in the condensate and the thermally excited quasiparticles. 
For a uniform Bose gas, our kinetic equation for the thermally excited quasi- 
particles reduces to that found by Eckern as well as Kirkpatrick and Dorfman. 



I. INTRODUCTION 



In a trapped, weakly-interacting Bose gas at T=0, the fraction of atoms that are excited 
out of the condensate is only a few percent [[J. As a result, the dynamics of the trapped Bose 
gas at low temperatures (compared to Tbec) is well described by the equation of motion for 
the macroscopic wavefunction $(r, t). This is the time-dependent Gross-Pitaevskii equation 
(GP) i 



d$(T,t) 
i- 



-^-V 2 r + U ext {r)+gn c (r,t) 
2m 



*(T,t), (1) 



dt 

where n c (r, t) = |$(r, t)\ 2 is the non-equilibrium density of the atoms in the condensate and 
U ext (r) is a harmonic trap potential (In this paper, we set h = 1). For a discussion of the 
properties of a dilute Bose gas at very low temperatures, only the s-wave component of the 
two-body interaction v (r — r') is important. Thus one can use the pseudopotential 

v(r) = gS(r), g = Aira/m, (2) 

where a is the s-wave scattering length of the true potential. The GP equation describes the 
motion of the condensate moving in the dynamic Hartree mean-field produced by the other 
atoms in the condensate and gives a closed equation for the order parameter $(r, t). The 
GP equation ([]]) provides a very accurate description of the static and dynamic properties of 
a trapped Bose gas at low temperatures T < 0AT BEC , as confirmed by many experiments 
in the last few years Q. In superfluid 4 He, the non-condensate fraction at T=0 is close 
to 90% 0. Thus in superfluid 4 He, one always has to deal with both the condensate and 
non-condensate atoms. Clearly a closed GP equation for $(r, t) like (|l|) is never valid in 
superfluid 4 He. 

At finite temperatures (say T > 0.5 Tbec)-, however, the number of atoms thermally 
excited out of the condensate becomes significant and the GP equation (|1|) is no longer 
sufficient. The simplest way to include the effect of the excited atoms on the condensate 
is to add the additional Hartree-Fock mean field Vhf = 2g ( n(r, t) produced by the non- 
condensate atoms (here n is the local non-condensate density). One immediately sees this 



new GP equation is no longer closed since it depends on the dynamics of the non-condensate 
atoms. 

To find the time- dependent non-condensate density, Zaremba, Nikuni and Griffin (ZNG) 
|| have used a quantum Boltzmann equation for the single-particle distribution function of 
the non-condensate atoms /(p, r, i) 

(3) 

coll 

Here, the thermally excited atoms are assumed to be well described by the single-particle 
spectrum ^ — h U(r, t) , where 

U(r, t) = U ext (r) + 2g [n c (r, t) + n(r, t)] (4) 

includes the self-consistent Hartree-Fock dynamic mean field involving the total time- 
dependent local density n(r, t). The right-hand side of @ describes the effect of collisions be- 
tween atoms on the time evolution of the distribution function /(p, r, t). In Bose-condensed 
gases, this collision integral has two distinct contributions 

= C 12 [f] + C 22 [f]. (5) 

coll 

Here, C 22 denotes the part of the collision integral that describes two-body collisions between 
non-condensate atoms. Above T BEC , this is the only term present. In contrast, C\ 2 describes 
collisions involving non-condensate atoms and one condensate atom. The role of C\ 2 is 
crucial since it couples the condensate and non-condensate components. 

ZNG also derived a generalized Gross-Pitaevskii equation that includes the effect of 
the collisions between the atoms in the condensate and the thermal cloud. Recently, we 
obtained the equations of motion derived by ZNG at finite temperatures in a more elegant 
way using the well-known Kadanoff-Baym Green's functions formalism Q. We also note 
that several other groups have also recently discussed the finite temperature dynamics of a 
trapped Bose-condensed gas 0-3, each using a somewhat different formalism. 

The kinetic equation (|3|) is valid in the semiclassical limit only: it assumes that the 
thermal energy is much greater than the spacing between the trap SHO energy levels 

3 
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df(p,r,t) 
dt 




{ksT 3> huo, where uq is the harmonic well frequency) as well as the average interaction 
energy (ksT ^> gn). ZNG have given a detailed derivation || of ([3]) at finite temperatures 
for a trapped Bose gas using the approach of Kirkpatrick and Dorfman || , who considered a 
uniform Bose gas. However, the ZNG theory is not applicable to very low temperatures be- 
cause the thermal excitations on which it is based do not include the collective (or phonon) 



part of the Bogoliubov spectrum E p = J (p 2 /2m) 2 + gn c p 2 /m. To make this generalization, 
one has to formulate a kinetic theory in terms of such Bogoliubov quasiparticle excitations. 
In the paper, we use the Kadanoff-Baym (KB) non-equilibrium Green's functions method 
[10| to derive such a generalized kinetic equation for the thermally excited Bogoliubov quasi- 



particles. To do this, we work within the second-order Beliaev-Popov approximation ||11|| . 

Kadanoff and Baym first formulated the method of deriving a kinetic equation for a 
normal interacting system using non-equilibrium Green's functions [|1(J. Kane and Kadanoff 



(KK) [|12| , |T3l| generalized this method to deal with a Bose-condensed gas, with the specific 



goal of using the resulting kinetic equations to derive the two-fluid hydrodynamics equations 
of Landau WM. An excellent review of the nonequilibrium real-time Green's functions and 



the generalized kinetic equation for the normal systems (non Bose-condensed) can be found 



in the book by Zubarev, Morazov and Ropke WE 



It is important to emphasize that although our analysis involves the non-equilibrium 
generalization of the Beliaev second-order self-energy used for systems in thermal equilib- 
rium, our work is quite different from the recent papers discussing the poles of equilibrium 
Green's functions within the second-order Beliaev approximation |IT|,|16|,|17| . In a very ele- 



gant formulation, Giorgini |L(J has calculated the quasiparticle energy and damping at finite 
temperatures in a dilute Bose gas in the collisionless regime by linearizing the equations of 
motion for fluctuations using the first-order dynamic Hartree-Fock-Bogoliubov approxima- 



tion. This leads, as expected, to the same excitations spectrum found by Shi and Griffin flT 
who calculated directly the poles of the single-particle equilibrium Green's functions using 
the second-order Beliaev self-energy diagram contributions. 

In the present paper, we use the second-order Beliaev approximation to discuss the 



non- equilibrium dynamics of a trapped Bose-condensed gas at finite temperatures. We use 
the second-order Beliaev self-energies with the lower order Bogoliubov excitation spectrum, 
including off-diagonal single-particle propagators, but we ignore the anomalous correlation 
function m. This last assumption defines what we call the Bogoliubov-Popov approximation 
|Tj|. It is important to note that, in this paper, we are primarily interested in the damping 
effect arising from the collisions between atoms. We do not calculate the second-order 
corrections in g to the quasiparticle spectrum or to the condensate chemical potential that 
are associated with the real parts of the second-order Beliaev self-energies. One can show 
that the real part of the second-order self-energies enter into the kinetic equation mainly 
through renormalized quasiparticle energy (see Eq. (6.3.77) in Ref. Hl5|). This suggests a 
simple way of extending the kinetic equations we derive by using an improved quasiparticle 
spectrum. 

We derive the kinetic equation for the distribution function for the thermally excited 
quasiparticles, as well as a generalized equation for the Bose condensate order parameter. 
The kinetic equation we obtain is the same as the one derived for a uniform Bose gas 
by Eckern QHJ in 1984, and by by Kirkpatrick and Dorfman (KD) in 1985 0. The KD 
derivation was based on a direct extension of the traditional method used to derive kinetic 
equations for classical gases, which obscured much of the physics. Moreover, KD did not 
explicitly derive equations of motion for the condensate degree of freedom. In a non-Bose 



condensed uniform gas, similar kinetic equations are derived in Ref. |15j using the related 
Keldysh formalism. 

The kinetic equation for the quasiparticle excitations which we derive in this paper, cou- 
pled to a generalized GP equation, provides a platform for studying different non-equilibrium 
aspects of a dynamics of a trapped, Bose-condensed gas at all temperatures, both in the 
collisionless and hydrodynamic domains. Linearizing our equations of motion around static 
equilibrium, one could calculate the density response functions which would exhibit collective 
mode resonances with a spectrum which goes past the generating Beliaev approximation. 
The high temperature limit of these coupled equations |||| has been recently used to study 



the dynamics of the condensate formation and growth in an inhomogeneous, trapped Bose 
gas ISO]. By taking the moments of the kinetic equation, the two-fluid hydrodynamic equa- 



tions have been derived [R . 2 1 



The present paper is a natural generalization of our two earlier papers based on the 
KB formalism. In Ref. |22||, we have derived kinetic equations within the full Hartree-Fock- 



Bogoliubov approximation but ignored collisions. In Ref. [§J, we derived a kinetic equation 
including collisions, but which was based on a simple Hartree-Fock particle-like spectrum 
and hence was not valid at very low temperatures. 

In Section II, we review the general equations of motion for the non-equilibrium Green's 
functions describing the non-condensate atoms as well as the equation of motion for the 
macroscopic order parameter. In Section III, we transform these equations to a local rest 
frame of reference where the order parameter $(r, t) is real, i.e., to a frame where the 
local superfluid velocity is zero. This naturally introduces the superfluid velocity and the 
local chemical potential as the spatial and time derivatives, respectively, of the phase of the 



order parameter |T2|]. We then specialize our equations of motion for the non-equilibrium 
Green's functions for the case of slowly varying external perturbations. We use the key 
assumption that all correlation functions vary slowly as a function of center-of-mass space- 
time coordinates but are dominated by small values of the relative coordinates. Following the 
Kadanoff-Baym approach |10|,|12|,^3| , we derive a generalized quantum Boltzmann equation 
for the frequency-dependent quasiparticle distribution function /(p, u>] R, T). 

In Section IV, we use this generalized KB quantum Boltzmann equation to derive the 
kinetic equation for quasiparticles at low temperatures, with the collision integral describing 
collisions between quasiparticles in the Bogoliubov-Popov approximation. In Section V, we 
discuss a general form of the local quasiparticle distribution function and introduce a new 
quasiparticle chemical potential situations in which the condensate and thermal excitations 
are not in diffusive equilibrium. In Section VI, we derive, in a self-consistent manner, a 
generalized Gross-Pitaevskii equation. In Section VII, we verify that our coupled equations 
exhibit the Kohn mode corresponding to the harmonic oscillations of the center-of-mass of 



the equilibrium condensate and non-condensate density profiles. 



II. EQUATIONS OF MOTION FOR NON-EQUILIBRIUM GREEN'S 

FUNCTIONS 



For convenience, we first review the KB formalism JTD] already used in our earlier work 
. In terms of quantum field operators, the many-body Hamiltonian (K = H — HqN) 
describing interacting Bosons confined by an external harmonic potential U ext (r) is given 
by: 



K 



2m 



1 f 

+ - J dvdr'ip\r)il)\r')v (r - r')T/>( r M r ')- 



(6) 



We separate out the condensate part of the field operator in the usual fashion p3|j24 



^(r) = (V»(r)) t + ^(r), 



(7) 



where = and {ijj(r)) t = <3>(r, t) is the Bose macroscopic wavefunction. The non- 

condensate (or excited-atom component) field operators ^(r) and ip\r) satisfy the usual 
Bose commutation relations. 

In a Bose-condensed system, the finite value of $(r,t) leads to finite values of the off- 
diagonal (or anomalous) propagators (^(l)^^')) an d (^(1)^(1')) . These must be dealt 
with on an equal basis with the diagonal (or normal) propagators, and thus it is convenient 
to work with a single-particle 2x2 matrix Green's function defined by | fL2| , p5[ 



g(l,l';U) = -i 



1 (7^(1)^(1')) CWM1')> ^ 



V 



(8) 



<T^(1)^(1')> (T^(l)^(l')) 
Here, T represents the time-ordering operator and we use the usual KB abbreviated notation, 
1 = (r,t) and 1' = (r',t'). We define g K and <p by 



g(l,l';U)=g > (l,l';U) t 1 >t v 
= g<(l,l';U) h<t v 



(9) 



Using (0), the matrix propagator in @ splits into two parts 



17) = 0(1,1'; U) + h(l,l';U). 



(10) 



(11) 



Here g is identical to (|8]), except that it involves the non-condensate part of the field oper- 
ators. The condensate propagator is given by 

h[l, 1 ; U) = — % 

$*(1)$*(1') $*(1)$(1') 

with (^ f (r)) t = $*(r,i). 

A very useful and elegant way of generating the equations of motion for both g and $ is 
to use functional derivatives with respect to weak external fields |T0|j25|j26|] . 

1 



/ 



dr 1 d2V t (l)C/(l,2)^(2)+ dv 1 V + (l>w(l) +^(lk 



(12) 



Here £7(1, 2) is an external generating scalar field non-local in space and time. It represents 
a perturbation in which a particle is removed from the system at point 1 and added at 
2. The symmetry-breaking fields r] ext and rf ext describe particle creation and destruction 



25,27|. All higher-order Green's functions can be neatly expressed as functional derivatives 



of single-particle Green's functions with respect to such generating fields. 

Following the Kane-Kadanoff (KK) analysis [p~2| , p~0[| , the Dyson-Beliaev equations of mo- 
tion for the real-time non-condensate propagators 0(1,1') can be written in the following 
2x2 matrix form 



dl 



00^(1, 1)-± HF (1,1) f (1,1') 



< 



dir(i,i)0 (l s i') 



dlS>(l,l)a(l,l'), 



(13) 



and 



< 



dW (1,1) ^(l, 1') - X HF (l, 1') 



dla(l,l)E>(l,l') 



< 

, 7 ;> 



dig (i,I)r(I,i'). 



(14) 



Here a(l, 1') and f (1, 1') are defined by the matrix elements 



a a p(l,l') = ~g>p(l,l')-g< p (l,l') 

r a ^(i,i') = 2^(1,10-2^(1,1'). 



(15) 



The non- equilibrium single-particle spectral density a a p(l, 1') will play a crucial role in this 
paper. In the above equations and elsewhere, integration over dl means integration over the 
coordinates (i"i,ti) and a trace over the matrix index a%\ and 5(11') = 5(r — r')5(t — t'). 

The single-particle self-energy which is involved in ( |I"3| ) and (|14]) has already been split 
into two parts 



m,i') = £" F (i,i') + £ c (i,i'). 



(16) 



The first-order Hartree-Fock self-energies are given by 



/ 



2n(l), m( 



1)\ 



5(11'), 



(17) 



\m*(l), 2n(l) j 

and S c is the second-order "collisional" Beliaev self-energy. The total density is given by 
n(l) = ig[\(l, 1 + ) + |$(1)| 2 = n(l) + n c (l) and the total anomalous density by m(l) = 
z<7i 2 (l, 1) + [$(1)] 2 - In addition, we define (see (|)) 



E c (l,l') =S>(1,1') t x >t v 

= E<(1,1') 



(18) 



In (|T3| ) and (Q), the inverse of the non-interacting 2x2 matrix Bose gas propagator go(l? 1') 
is defined by 



iT 3-^T + 7. UeBt^Tx) + Ho 

oti 2m 



5(1,1'). 



(19) 



We note that the equations in (|T3| ) and (|T4| ) already have the "structure" of a kinetic 
equation such as (|). The Hartree-Fock part of the self-energy has been included into the 
left-hand side of (O) and (|HD , giving the mean- field contribution to the "streaming" term. 



The second-order self-energy describing binary collisions is included on the right-hand side 
of ( |T3"D and (|i4|), and it will be shown to give rise to collision integrals in the quasiparticle 
kinetic equation we derive. 

In this paper, we work with the second-order self-energy S c as given by the Beliaev 
(gapless) approximation |TT],[rS],|25| . The advantage of the Beliaev approximation is that 
the non-condensate single-particle Green's function exhibits the correct quasiparticle spec- 
trum (phonon-like in the long- wavelength, uniform gas limit). In the second-order Beliaev 

> 

approximation, the second-order self-energy S < is given by {| 



£<(!,!') = -\v{U)v{2l')g<{ll' 



^>(23)^<(32) + g>(23)h(32) + h(23)g<(32) 
~g > (23)h(31') + /i(23)S<(31') + g>(23)g<(3l' 
h(ll')g>(23)g<(32) + 2h(12)g>(23)g<(3l' 



-v(13M21')0<(12) 
- -u(13)i;(2l') 



(20) 



The equation for the condensate can be written (see Ref. 0] for more details) in terms 
of the 2-component order parameter G\/2{1) = \/— where \& is defined as 



*(1) 



/ ^(1) ^ 



(21) 



One finds that equation of motion for the G\ji is given by 



J dlgo\l, 1)G 1/2 (1) = V=ir){l) + v^Wl), 



(22) 



where the so-called condensate source function 77 is defined by the three-field correlation 
function 



V /= I^(1) = -J d2V = iv(12)(TV(l)&(2)V(2)). 



(23) 



The external particle-source fields are defined in (|T2|), with 



VextO-) 



Vext 



Vtxti 1 ) 



(24) 
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The exact coupled equations of motion fljjD, (0) and (^) are the starting point of our 
analysis. The external generating fields U and r] ext will be left implicit in the rest of this 
paper. 

Using (|7|), one can also decompose the three-field correlation function involved in 77 
defined in (^). For example, one has 

(7^(1)^(2)^(2)) = $(1H(2) + $(2)<7V(1#(2)) 

+ $*(2)(7Y(1M2)> + (2^(1)^(2)^(2)). (25) 

In the first order Hartree-Fock-Bogoliubov (HFB) approximation, one neglects the three- 
field correlation function {^ijr ijj) for the non-condensate atoms. In this approximation, (j25l) 
only involves the condensate density n c and the two-field correlation functions n and rh, as 
defined earlier. In this paper, in contrast, we keep the non-condensate three-field correlations 
but will eventually set rh = (the Popov approximation p5| , p!8| ). 

Equation (p^) can be rewritten in terms of a condensate self-energy function 5, defined 

by 

J dlS(l,l)h(l,l') = V^v(l)G\ /2 (l'), 
with the condensate propagator h given in ([TTJ) . In place of (|2"2"D, we have 



(26) 



/ dl feo-^l, 1) - S HF (1, 1)1 h(l, 1') = f dl(S>(l, 1) - 5<(1, l))/i(l, 1') 



(27) 



where, as before, the mean- field contributions (see fl25|) ) are included in the Hartree-Fock 
part of the condensate self-energy S HF 

-)+n(l), m(l) ^ 

5(11'). (28) 

^m*(l), n(l)+n(l) y 



/ 



S HF (11') 



9 



The functions S < on the right-hand side of (|27| ) contain the second-order contributions to 
the condensate self-energy H 



£<(1,1') = -^(13M21')r(ll') 

> 



-w(13>(21')^ < (12 



£>(23)£<(32) 
r(23)r(31' 



(29) 



11 



As we will show in Section VI, these contributions give rise to a dissipative term in a 
generalized GP equation. If we recall the definition for the condensate propagator h in ( |TT|) 
and use the explicit form for the condensate Hartree-Fock self energy S HF in (|28|), we obtain 
a generalized GP equation for the macroscopic wave function $(r, t), namely 

-J- + H - tWri) +l*o-9 Ml) + 2n(l)) $(1) = £m(l)$*(l) 
+ /* dl (S> - S<) (1, T)$(l) + f dl (S> - S<) (1, (30) 

J —oo J -co 

Initially, the equation of motion for the macroscopic order parameter fl2"2"|) was given in terms 
of the condensate three-field source function 77 in (|23|) . We have rewritten the equation for $ 
in terms of the condensate self-energies S defined by (p6| ) because, as we will see in Section 
VI, S is more convenient to work with than the condensate source function rj. 

III. GENERALIZED KINETIC EQUATION 

In general, the order parameter $(r, t) in a Bose fluid is complex. It is often written in 
terms of the condensate amplitude and phase 



$(r,t) = ^ c (r,t)e 4 ^. (31) 

To derive a generalized kinetic equation, the rapid oscillations of the phase in ( |3~TD cause 
difficulties. Following KK, we first gauge transform (^3|) and ([14]) to the local rest frame 

in which the superfluid velocity is zero. This corresponds to removing the phase of the 

< 

macroscopic wavefunction. The required gauge transformations on h(l,V) and g (1, 1') are 



12 



h'(l,l') = e-^ T<3) h(l,l')e^ T(3) ., 



^(1, 1') = e~^ (3) f (1, l')e^ (3) , (32) 

where r( 3 ) is the Pauli spin matrix. The physical interpretation of (^) is that it involves 
a transformation to a coordinate system in which non-condensate atoms are moving with 

12 



average velocity v s relative to a stationary condensate. For example, the transformation 



(|32|) gives 



(l,l') = e-W)-^ 1 <(l,l'), 



^(1,10 



-iyJn c (l)n c (V). 



(33) 



One sees that, in the local rest frame (denoted by a prime), the order parameter $ (r, t) is 
real (see (0)). 

Equations ( |13|) and fll4j) remain unchanged in form after this transformation to the local 
rest frame as long as g^ 1 is replaced by (compare with ( |19D in the lab frame): 

d 86(1) 1 



5(1,1')- (34) 



We recall that the superfluid velocity v a (R, T) and the local chemical potential /i c (R, T) are 
defined as the spatial and time derivatives of the phase fll^] , namely 



mv s (R,T) = V R #(R,T) 
<90(R,T) _ 



<9T 



1 



/x c (R, T) - /i + ^mf s 2 (R, T) 



(35) 



where, in the lab frame, the condensate wavefunction is given by $(R, T) 



n (R, T)e ie ( n ' T \ Therefore, the gauge transformation changes the momentum p — > 
p — mv s as expected for the momentum in the local rest frame (when we Fourier trans- 
form, the gradient operator in (|34|) becomes the momentum). 

In Section II, we have written down the equations of motion (pT3|) and ( |14D for the non- 
equilibrium real-time Green's functions. We now want to use these to derive a kinetic 
equation for the quasiparticle distribution function. We recall that (|13|) and (|14]) are matrix 
equations. In addition, (O) involves differential operators with the respect to the coordinates 



(r, t) and (|T^) involves derivatives with respect to the coordinates (r',t f ). Since our single- 
particle Green's functions are functions of both coordinates (1,1') , one has to find a way to 
combine both equations to derive a single kinetic equation for a quasiparticle distribution 
function. We will discuss later how these Green's functions are related to the quasiparticle 



13 



distribution function. However, we emphasize that we will need to use both ( |i~3"D and fll4"D 
to derive a kinetic equation for quasiparticles. 

In the Kadanoff-Baym procedure P,|i~0|| one rewrites the equations of motion in terms of 
relative and center-of-mass space-time coordinates, defined by 

+ + + T> r l + r l' rri h+tll . . 

r = ri — iv, t = t\—ti>\ R = , 1 = . {So) 



In the literature, this is sometimes known as the Wigner representation |T5|J28[] . It allows 
one to separate variables describing "slow" and "fast" processes in the system. In thermal 
equilibrium, the Green's functions are only functions of the relative space-time coordinates 
r and t, and moreover are sharply peaked about r = and t — |10| . Since we assume that 
the external disturbances are slowly varying in space and time (with a wavelength much 

larger than the thermal deBroglie wavelength), we expect that these slowly varying external 

> 

disturbances will not change this dependence of ^ < (r, t;R, T) on small values of r and t. 
Therefore, our non-equilibrium correlation functions (like g, E, etc.) are assumed to be 
dominated by the small values of relative coordinates (r, t) (equivalently, by high momenta 
and frequencies in the Fourier transforms), but vary slowly as functions of the center-of-mass 
coordinates (R, T). 

Using these key properties of the non-equilibrium correlation functions to simplify the 
equations, we now write ( |I~3| ) and (14]) in terms of the center-of-mass and relative coordinates 



T^JTBIP^ ], and take the trace of the resulting matrix equation to obtain: 

C u g< (r, t; R, T) + t 22 ~g< (r, t; R, T) — g 
+ J°° dfdtTr\E > {r-f,t-f)g~ < {f,^-E < {r-f,t-^g>(f, : t)^. (37) 



((r • V R + t-^)m)g< + (I r • V R + 1-^ J m* J gf 2 



As usual, for simplicity, the (R, T) dependence of the g, £ and m is left implicit. The 
operators L\\ and £22 on the left-hand side of (|37|) are defined as 



d d 1 

Ai = + (r • V R + t— ) (/i(R, T) - U eff {R, T)) + -V R ■ Vr 

+ % [ (r ■ V R + ^)v s (R, T) ) ■ V r + zv 8 (R, T) ■ Vr + *V R • v s (R, T) 



dT' 
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d d X 

£22 = -igf + (r ■ V R + t— ) (MR, T) - U eff (R, T)) + -V R • Vr 

-i|(r-V R + ^)v s (R, T)^ ■ V r - iv s (R, T) ■ V R - zV R ■ v s (R, T), (38) 
where the effective dynamic HF field f/ e //(R, T) is given by 

C/ e// (R, T) = U ext (R) + 2<?(n c (R, T) + n(R, T)). (39) 
We emphasize that in the expansion for the small values of relative coordinates (r, t), we 



did not keep all terms of order d/dT and Vr, in (p7|). These additional terms that we have 
neglected contribute to the many-body renormalization effects, i.e. how the two-particle 
interaction changes the dispersion relation of the quasiparticles due to terms of second 
order in g. Such corrections involve the real part of the second-order Beliaev self-energies. 
The Bogoliubov-Popov quasiparticle approximation we use for the spectral densities a a p 
in fllS!) do not include such second-order effects. In the present paper, we concentrate on 
the damping effects associated with the collisional self-energies S < on the right-hand side of 
fl37p. For further discussion of the KB formalism related to going past the simple Bogoliubov 



quasiparticle approximation, see Ch. 9 of Ref. JTO | and Ch.6 in Ref. | 15[ 



The double Fourier transform of (|37|) gives 



dfn* dg^~ dirt dg^ 

Cn9ii + £22^2 - #V R m* ■ Vp<?f 2 - #V R m ■ Vp^ + 9^r^- + 9 gT g ^ 

= Tr (±>(p, 00; R, T)g<(p, uj; R, T) - S< (p, u; R, T)g>(p, w; R, T)) , (40) 



with 9ap = 9ap(P> R, T) and 

d _ d _ d 

Cn = — + V p [e p + p ■ v s ] • V R - V R [e p + p ■ vj • V p + — [e p + p • v fl ] — 

d d d 

£22 = --^ + V p [i p - p • v s ] • V R - V R [t p - p ■ v s ] • V p + — [i p - p ■ v s ] — . (41) 

The result in (^DJ) gives an equation closely related to the quasiparticle kinetic equation we 
are trying to derive. Here, e p is defined by 

2 

e p (R, T) = 2- + U ext (R) + 2gn(R, T) - /i c (R, T). (42) 

15 



In Section VI, we shall see that the condensate chemical potential // C (R, T) is given by 



He = V , + U ext (R) + U(R, T)+g [2h(R, T) + n c (R, T)\ . (43) 

2mJn c (R,T) 



A kinetic equation for thermally excited atoms in a trapped Bose gas can be written 
in terms of distribution functions for either atoms or for quasiparticle excitations. In our 
earlier papers |4.22], we transformed the equations of motion for a non-equilibrium Green's 



functions at high temperatures into a kinetic equation for a single-particle distribution func- 
tion /(p, R, T) describing the non-condensate atoms. The latter are assumed to have a 
Hartree-Fock spectrum. If one wants to use a more realistic spectrum valid at low temper- 
atures, it is much more convenient to work within a quasiparticle picture. In the theory of 
Bose-condensed trapped gases, one introduces quasiparticles by expressing the quantum field 
operators for the non-condensate atoms as a coherent superposition of creation and annihi- 
lation operators for Bose quasiparticles, with the weights given by the usual Bose-coherence 
factors u and v 

^(R, T) = J2 [u i (R)a i e- iEiT ^ + v*(R)a\e lE * T/n ] . (44) 

i 

In the semiclassical approximation, ([14]) becomes 

^(R, T) = J-^y 3 [u p (R)a p e-^ h + V ;(R)6^ T ^] . (45) 

Here a p and a p are the Bogoliubov quasiparticle creation and annihilation operators, re- 
spectively, which obey the usual Bose commutation relations. One can see that creating an 
atom with momentum p is equivalent to creating a quasiparticle with momentum p with 
amplitude u p and at the same time, destroying a quasiparticle with momentum — p and 
amplitude v p . The quasiparticle distribution function is given by the statistical average of 
the quasiparticle operators, i.e., /(p) = {aJ p a. p ). We recall PJ2^| that distribution function 
for atoms f at (p,uj] R, T) is directly related to the diagonal Green's function, namely 

f at (p, u; R, T) = -igf x (p, u\ R, T). (46) 
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In terms of quantum field operators, is given by (see Ch. 9 of Ref. fLOf ) 

~g< (p, W ; R, T) = i J drdte-^ M (ft(R - T - ^(R + T + t)). (47) 

The usual Boltzmann equation is expressed in terms of the Wigner distribution function 
fw(p, Rj T) j|. This limits the description to the semiclassical approximation because it is 
assumed that the position and momentum of the particles can be defined simultaneously. 
In order to use this kind of distribution function for quantum systems, it is necessary to 
perform some type of averaging in order to remove effects due to the uncertainty principle. 
In our work, we want to derive a kinetic equation for the quasiparticles which is valid at all 
temperatures and therefore the semiclassical approximation will no longer be valid. To in- 
clude the quantum effects, we introduce the quasiparticle distribution function /(p, u>; R, T) 



with an additional variable to in the following way ||12| , p!5| , p8| , p9 
QapiP' w ! R ' T ) = ia ap(p, w; R, T)/(p, u; R, T) 

^(p, u; R, T) = ta af3 (p, lu; R, T) [1 + /(p, lu; R, T)], (48) 

where the spectral density a a/ 3 is defined in (^5]). Using the Bogoliubov- Popov approximation 
for the spectral density in (|54D , one can use ( ^8|) to obtain the well-known relation between 
the quasiparticle distribution function /(p, u;~R, T) and the atom distribution function f at 
in (|46|), namely 

f at (p, u- R, T) = («J(R, T) + Up(R, T)) /(p, R, T) + i; p 2 (R, T). (49) 

The semiclassical Wigner distribution function is obtained by taking the frequency integral 
of the atom distribution function defined in ([49]). 

One can show, using ([13]) and (0), that the spectral function a a p in ( [15] ) satisfies 

C n a n + t 22 a 2 2 - g^nm* ■ V p a 12 - gV R m ■ V p a 21 - 9^r^^ ~ = °- ( 50 ) 



Using ( ^8|) and (|50|), one can rewrite the kinetic equation (^) for ^< to obtain a new kinetic 
equation specifically for the quasiparticle distribution function /(p, o>; R, T) in the following 
form: 
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+ a ^ d -w% + fl2 w£ = /Tr(s>a) - (1 + ^ Tr ( s<a )- ( 51 ) 

We notice that fl5T] ) includes terms involving 9/ 9a;. The additional variable a; in /(p, u; R, T) 
results in new streaming terms on the left side of flSTp. These terms are not present in the 
semiclassical kinetic equation (which is obtained from (pT|) integrating over a;), which shows 



that the terms involving d/du are of quantum origin |15| , |28| , |29 . 

Eq. fl5l]) is the most general form for a kinetic equation for the quasiparticle distribution 
function / within our model. To derive (pi]), we have only assumed that the external 
disturbances vary slowly in space and time, and therefore all relevant physical quantities 
vary slowly as function of center-of-mass coordinates (R, T) defined in (0). The other 
assumption that we made is that one can introduce a quasiparticle distribution function 
/ through the definition in (f£8"D- Of course, at this stage, one could say that we are only 
replacing one unknown function with another. The generalized kinetic equation for a Bose- 
condensed system ( |5"T| ) was first derived by Kane |Tj|. From (|5TD, we see that the general 



structure of the collision integral I has the following form [p^ ,|l3|] 

7[/(p, R, T)]= fp- \fTr{±>a) - (1 + f)Tr(E<a)} . (52) 

> 

Using (^) in the the general expression for the non-equilibrium Beliaev self-energies S < 
given in ( p0|) , one can prove that the collision integral given by (|52|) conserves momentum 
(see Appendix for details), 

/dppI[/(p,R,T)]=0. (53) 

One can also prove that ( p^ ) conserves energy as well. To prove this, however, we need to 
work within a specific approximation for the single-particle spectral density a aj a(p, u; R, T). 

IV. KINETIC EQUATION IN THE BOGOLIUBOV-POPOV APPROXIMATION 



In this section, we use the results of Section III to derive a quasiparticle kinetic equation 
within the Bogoliubov-Popov approximation. More precisely, this means that we will use 
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the spectral densities with the Bogoliubov-Popov quasiparticle excitation energies |T0,|T3|J23 



a n (p,u;R,T) = 2ir[u 2 p 5(uJ - v s • p - E p ) - v 2 p 5(u - v s ■ p + E p )} 
a 12 (p,a>; R,T) = -2ttu p v p [5(uj - v s • p - E p ) - 5(u - v s • p + E p )\ 
a 2 i(p, u; R, T) = ai 2 (p, u;; R, T) 

a 22 (p,a;;R,T) = — on(— p, — a>;R,T). (54) 
Here, the Bose- coherence factors w(R, T) and t>(R, T) are given by ]23| 



and the quasiparticle energy E p is given by 



E P (R, T) = Je?(R, T) - (^ C (R, T)) 2 . (56) 



We emphasize that the spectral densities in (|54|) could be derived in the quasiparticle ap- 
proximation from the general equations of motion for the Green's functions as it has been 



shown in [T3[]. We simply start with them as input into our general formalism. 

In the Thomas- Fermi approximation one neglects the quantum pressure term in (f43D, 
in which case the quasiparticle energy E p reduces to the usual Bogoliubov excitation energy 



E P (R, T) = + 2gn c (R } T)e p , (57) 

where e p = p 2 /2m. We note that spectral densities a a p(p, uj; R, T) in ([54]) exhibit both 
positive and negative energy poles. In the Hartree-Fock approximation used in our earlier 
work [Q, u p = 1 and v 2 = 0. Physically, ( |54"D corresponds to the assumption that the thermal 
cloud can be considered as a gas of weakly-interacting single-particle excitations with the 
excitation energy given by (0). One can check explicitly that ( |54]) do satisfy the general 
equation of motion given in (0). Note that in the literature, v p is sometimes defined with 
the opposite sign, such that u p v p in ([55]) is negative. 

Substituting the spectral densities fl5"4] ) into Q5I|), we can now derive a kinetic equation 
for the quasiparticles. After lengthy algebra, we obtain 
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00 du 
-oo 2tt 



5(u — v 8 -p - E p 



2i 

dT 



+ V p (E p + v s ■ p) • V R / - V R (£ p + v s • p) • V p / 



<9 <9f 

^ (25 P + v s • p) •£ + ^ ((1 + /)£< - /S>) + „J ((1 + /)£< - /S> ) 



W ((1 + /)(S< +S< 



21/ 



—5(uj - v s • p + Ep 
-oo Air 



/(£> + £>))] 

^ + V p (-E p + v s • p) • V R / - V R (-E p + v s • p) • V p / 



+ w (-Ep + v s • p) -£ + vl ((1 + /)£< - /E>) + u\ ((1 + /)E< - /£> ; 
- W ((1 + /)(£< + £<) - /(E> + E>))] = 0, 



(58) 



where Eq^ = E a/ 3(p,u;; R, T). This is a kinetic equation for the frequency dependent quasi- 
particle distribution function / expressed in terms of an integral over both positive and 
negative energy poles. If we recall the expression for Bose coherence factors u and v given 
by (0), we note that the second term in fl5"S| ) is the same as the first term in fl5"8]) if we 
replace — E p with E p . Therefore, it is sufficient only to consider the first term to obtain the 
kinetic equation for the quasiparticle distribution function defined by 



f qp (p, R, T) = f(p, u - v s • p = E p - R, T). 



(59) 



We obtain finally 



df 

+ V p (E p + v s • p) • V R / W - V R (E p + v s • p) • V p f qp ] = I[f qp }. (60) 



Here the collision integral I[f qp (p, R, T)], defined in (|52|) , becomes 



I[/ ?p (p,R,T)] 



2tt 



«j ((1 + f qp )E< - / W E>) + d ((1 + /„,)£< - f qp T> 2 ) 



+ UpVp ((1 + / ?P )(E< + £<) - / ?P (E> + E> ))] 



(61) 



To evaluate the collision integral in ([H]), we need to choose a specific approximation 
for the second-order self-energy Y> a p. Here, we use the second order Beliaev approximation 
given by (120), the Fourier transform of which is 



£<(p,^;R,T) = -^ 2 



dpiduj t 
(2vr) 8 



S(oj + loi 0J2 o; 3 )5(p + pi - P2 - P3) 
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> 

9 (P2,w 2 )7V 



< 



> 



g (pi,uJi)g (p 3 ,^3 



> 



< 



> 



+ 2# (p 2 ,^2)5 f (Pl,Wi)p (p 3 ,^3 



> 



+ 9 (p2,u 2 )Tr 

+ h(p2,uj 2 )Tr 

> 

+ 2g (P2,W2) 



> 



< 



h(p 1 ,u 1 )g (p 3 ,u 3 )+g (pi,Ui)h(p 3 ,u 3 



< 



> 



9 (Pi,vi)g (p 3 ,^3 



< 



> 



+ 2h(p 2l uj2)g (pi,wi)g (P3,^3 



> 



< 



h(p 1 ,u)i)g^{p 3 ,u 3 )+f{p 1 ,u 1 )h{p 3 ,u 3 ) . (62) 

As usual, the (R, T) dependence of the functions T,,g and h on the right-hand side has 
been suppressed for simplicity of notation. The quasiparticle energy E P (H,T) in (|60[) is 
the energy of the quasiparticles in the local rest frame (v s = 0). The Beliaev second-order 
expression Q52] ) consists of two kinds of contributions: (1) Terms that include both the 
condensate propagator h and the non-condensate propagators g; (2) Terms that include 
the non-condensate propagators g only. The first type of contribution will give rise to the 
collision integral that describe collisions that include one condensate atom interacting with 
the thermally excited quasiparticles. As in earlier work |||],|[|, we denote this part of the 
collision integral as C\ 2l indicating that we go from 1 thermally excited quasiparticle (and one 
condensate atom) to 2 thermally excited quasiparticles. The second type of contribution in 
(|62|) only includes non-condensate propagators. We denote this part of the collision integral 
as C22, indicating that it describes collisions where 2 thermally excited quasiparticles are 
scattered into 2 excited quasiparticles. At low temperatures, when the number of thermally 
excited quasiparticles is small, we can neglect the C 22 collision integral relative to G\ 2 . 

From (|53|), it follows that both C\ 2 and C 22 conserve momentum (see Appendix for 
details), i.e., 



dppCi 2 = 
dppC 22 = 0. 



(63) 



In addition, using fl54|) in ( |62]) one can show that the collision integral in ([52]) conserves 
quasiparticle energy E p and therefore both C\ 2 and C 22 will satisfy the conditions 



/ dpE p C l2 = 



21 



I dpE t 



pC 22 



(64) 



One can show that for slowly varying external disturbances, the condensate propagator 
in (|TTD can be approximated by 



h(p, lu; R, T) = n c (R, T)(2ir) 4 5(p)5(u) 



(65) 



V 1 x / 



To evaluate the collision integral C\ 2 in terms of the Bose coherence factors u and v, we 
need to substitute ([54]) and (|65|) into d62|). One can simplify (|62"D greatly using the following 
exact symmetry relations; 



< > 



5f 22 (p, w; R, T) = ^ii(-p, -w; R, T), 

< > 

sfe(p, w; R, T) = <?f 2 (-p, -w; R, T), 

< > 
fifa(p,w;R,T) = <7n(-p, -w;R,T). 



(66) 



After some algebra, one finds the following expressions for the non-equilibrium self-energy 
(considering now only the C\ 2 terms, which include one condensate propagator): 



S<( P ,.;R,T) = V/^?- c (R,T) 



[P2, uJ 2 )9n(p ~ P2, w - w 2 ) 



+ 4^n(p2, uJ2)9ii(V2 - p, w 2 - w) + 8gf 2 (p 2 , £*> 2 )<7n(p - p 2 , a; - ^ 2 



+ 4<? 1 < 2 (p 2 , W2)^i < 2 (p - P2, a; - w 2 ) 



(67) 



£<(p,w;R,T) = -# 2 ^ d ^ 2 ^ 2 ra c (R,T) 



> > 

4^1 < 2 (P2, UJ 2 )gf 2 (p - P2, W - W 2 ) 



+ 4^h(p 2) uj 2 )gn(p 2 - p, w 2 - w) + 8gf 2 (p 2 , u 2 )g^ l {p 2 -p,u 2 - to) 



+ 2fi' 1 > 1 (- P2 , -uj 2 )g^(p 2 -p,u 2 -u 



(68) 



£<(p,u,;R,T) = -, 2 /^?n c (R,T) 



(2tt)< 



> 



> 

1 < 2 (P2, ^2)^i < 2 ( P - P2, W - W 2 ) 



> > > < 

+ 4^ 1 < 2 (p 2 ,cj 2 )^ 1 < 1 ( P -p 2 ,u;-^ 2 ) +4^ 1 (p 2 ,a; 2 )^ 1 > 1 (p 2 -p,w 2 - 

> < 
+ 4<?f 2 (P2, w 2 )^n(p 2 - p, w 2 - w) 



(69) 
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£<(p,^R,T) = £<(p,^R,T). 



(70) 



Again, the (R, T) dependence of the g's is suppressed on the right-hand side. One notices 
that these expressions have the same structure as the thermal equilibrium ones obtained by 



Shi and Griffin, as well as others |11]J16| , ^7|| . This is expected since the whole structure of 
our theory is only valid for a systems slightly perturbed from thermal equilibrium, with all 
physical quantities assumed to vary slowly as functions of center-of-mass coordinates (R,T). 
The entire KB formalism reduces to the usual equilibrium self-energies in the appropriate 
limit, which is one of its strengths. 

One can show using the general properties of the non-equilibrium Green's functions 
in fl66|) in conjunction with (|48|) that the quasiparticle distribution function /(p, u>; R, T) 
satisfies the exact relation 



/(-p, -w; R, T) = —(1 + /(p, w; R, T)). 



If we introduce following standard abbreviations for the Bose-coherence factors 



Ap — itp, Bp — Up, Cp 



-UpVp, 



the self-energies in (|6T|)-(|69"|) can be written as ( using ([H|)) 

dp 2 du 2 



£<( P ,u;;R,T) = //^n c (R,T) 



/ :i+/i)(i+/ 2 )^ 



v 



/i/ 2 



[{2A X A 2 + 8A X C 2 + 4CiC 2 + AB X A 2 ) 5{u 2 - E 2 )5(u x - E x ) 

- (2B X A 2 + 8B X C 2 + 4CiC 2 + 4A X A 2 ) 5{uo 2 - E 2 )5{oo x + Ex) 

- (2A X B 2 + 8A X C 2 + AdC 2 + AB X B 2 ) 5(u 2 + E 2 )5(u x - E x ) 
+ {2B X B 2 + 8B 1 C 2 + AC X C 2 + AA X B 2 ) 5{u 2 + E 2 )8{u x + E x )\ 



E%(p,u;R,T)=g 2 J 



dp 2 du 2 

—n c (R, T) 



( {l + h){l + f 2 )\ 



{2<Kf 



\ 



J 



[(6dC 2 + AA X C 2 + AB X A 2 + AB X C 2 ) 5{uj 2 - E 2 )5(u x - E x ) 



(71) 



(72) 



(73) 
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- {AB X C 2 + AA X C 2 + 6C X C 2 + AA X A 2 ) 5{u 2 - E 2 )5(co x + E x ) 

- (4AiC 2 + 4BiC 2 + 6CiC 2 + 45x5 2 ) 5(^ 2 + £ 2 )%! - E x ) 
+ (4Aifl 2 + 4 J B 1 C 2 + 6C X C 2 + 4A 1 C 2 ) 5(o; 2 + E 2 )8{u x + £1)] 



(74) 



and 



£< (p, cj; R, T) = g 2 J -^^ n c (R, T) 



/ (l + /i)(l + / 2 )^ 



V 



/i/a 



/ 



[(25^2 + 8£iC 2 + 4CiC 2 + AB X A 2 ) 5{u 2 - E 2 )6{u 1 - E x ) 

- (2A X B 2 + 8A X C 2 + AC X C 2 + AA X A 2 ) 5(u 2 - E 2 )5(u x + E x ) 

- (2B X A 2 + 8B X C 2 + AC X C 2 + AB X B 2 ) 5(uj 2 + E 2 )5{u x - E x ) 

+ (2A X A 2 + SA X C 2 + AC X C 2 + AA X B 2 ) 5(u 2 + E 2 )6{u x + E x )\ , (75) 

where p — p 2 = pi and u — u 2 = u x . 

Using these results, we can finally evaluate the C X2 collision integral given in (|6T|) 

C X2 [f] = 2g 2 n c (R, T) J [(1 + f)f x f 2 - f(l + f x )(l + f 2 )) 8(p - Pl - p 2 ) 

((«! - v x ){u p u 2 + v p v 2 ) + (u 2 - v 2 )(u p u x + v p v x ) - (u p - v p )(u x v 2 + v x u 2 )f 5(E p - E x - E 2 ) 
+ 2 {{u x - v x )(u p u 2 + v p v 2 ) + (u p - v p )(u x u 2 + v x v 2 ) - (u 2 - v 2 )(u p v x + u x v p )) 2 5(E P + E x - E 2 ) 
+ {v p u x u 2 + u p u 2 v x + u x v 2 u p - v p v 2 u x - v x u 2 v p - u p v 2 v x ) 2 5(E p + E x + E 2 ) . (76) 

The last term in QTB| ) clearly vanishes because of the energy delta function. We recall that 
all u's, u's and the quasiparticle energy E p in ( |76|) have an implicit (R, T) dependence. If 
we change pi — > — pi in the second term in (|76|), and use (|7l"l), we can simplify (|76|) slightly 
to obtain 



C X2 [f] = 2g 2 n c (R,T) J 



^%[(l + /)/i/ 2 -/(l + /i)(l + / 2 )] 



(2*) s 



((u x - v x )(u p u 2 + v p v 2 ) + (w 2 - v 2 )(upU X + Wp^i) - (it p - v p )(u x v 2 + ^iM 2 )) 2 
x<5(p - pi - p 2 )S(E p - E x - E 2 ) 
- 2 ((u x - v x )(u p u 2 + VpV 2 ) + (u p - v p )(u x u 2 + v x v 2 ) - (u 2 - v 2 )(u p v x + u x v p )) 2 

x5(p + Pl - P2 )5(Ep + E X - E 2 )] . (77) 
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The first term in (|77j) describes the decay of an excitation with momentum p into two exci- 
tations with momenta pi and p 2 . At T = 0, this is the only scattering process possible since 
there are no thermal excitations. The second term describes an excitation of momentum 
p absorbing a thermal excitation of momentum pi, leaving an excitation with momentum 
P2 = P + Pi- This form of the collision integral was first written down by Eckern ||19||, and 



shortly after Kirkpatrick and Dorfman || gave a more detailed derivation. Here, we have 
used the Kadanoff-Baym approach to give a cleaner derivation of C\ 2 , in a form which is 
also valid for a trapped Bose-condensed gas. 

After some algebra, one can also rewrite (|77|) in the following more compact form || 



C 12 [/] = 2g 2 n c (R, T) J dpi ff Pa | A(2, 3; 1) | 2 5( Pl - p 2 - p 3 )5(E 1 - E 2 - E 3 



dpidp 2 dp 3 

[ S (p - Pl ) - 5(p - pa) - 5(p - p 3 )] [(1 + f x )hh - /i(l + / 2 )(1 + h)] ■ (78) 



Here the scattering amplitude in \A\ 2 is given in term of the Bose coherence factors u and v 
A(2, 3; 1) = (u 3 - u 3 )(itiit 2 + v x v 2 ) + {u 2 - ^ 2 )(^iM 3 + ^1^3) - («i - fi)(w 2 f3 + v 2 u 3 ). (79) 

The first term in (^) is equivalent to the first term in (|77|), while the other term in ([77| ) is 
equivalent to the second and third terms in (f78T). 

In conclusion, the kinetic equation we have derived for thermally excited quasiparticles 
is given by 



d 

— + V p {E p + v s • p) • V R - V R (E p + v s • p) • V p 



f qp (p,R,T)=C 12 [f qp ] (80) 

with C\ 2 given explicitly by (|77|), or equivalently, (|78|). The derivation of this equation is 
the main result of this paper. 

To remove the rapidly varying phase of the order parameter real, we have gauge trans- 
formed to the local rest frame where the condensate is at rest. Hence, the energy of the 
thermally excited quasiparticles is measured relative to this local frame. Since the thermal 
excitations are moving with the superfluid velocity v s relative to the condensate, the energy 
of quasiparticles measured relative to the condensate is E p + v s -p. Therefore, the expression 
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Ep + v s • p in the streaming term on the left-hand side of (|80"D is expected. Similarly, if we 
denote the quasiparticle distribution in the coordinate system where the quasiparticles are 
at rest by /(p, u; R, T), the quasiparticles moving with the velocity v s relative to the con- 
densate will be described by the distribution function /(p, u — v s ■ p; R, T) = f qp (p, R, T), 
PaP- 



as occurs m 



If we use the frame of reference where the quasiparticles are at rest, the streaming term 
will include the energy of the quasiparticles only ( i.e., the v s • p term in (^) will not be 
present). This lab frame of reference is used in the work of Zaremba, Nikuni and Griffin ||. 



To understand the C12 collision integral in (|77|) better and the corresponding scattering 
processes that it describes, it is useful to consider a few limiting cases for a uniform gas. We 
define p\ = 2mgn c as the characteristic momentum for the crossover between the linear and 
the quadratic part of the quasiparticle spectrum (p = hk = where £ is the healing 

length). We then consider the following special cases: 

1) If all momenta Pi 3> po, then the quasiparticle spectrum E p defined by ([57]) is equal to 
a single-particle spectrum e p given in (pE2"p. Moreover in this limit, it follows from (|55|) that 
u — > 1 and v — > 0. Hence, the scattering amplitude A in ( |T5D becomes unity and the collision 
integral in ( |75| ) then reduces to the one recently derived by Zaremba, Nikuni and Griffin m 
using a different approach. Clearly, this approximation is only valid at finite temperatures 
where the dominant excitation spectrum is described by the Hartree-Fock single-particle 
spectrum in (|42|). 

2) In the opposite limit, when all three momenta Pi are small, one can expand the Bose 



coherence factors u and v in the following way [30,31 



™ X 1 / 2 1 / w \ 1/2 
gn c \ i 1 / h, v \ 



Up ~ I I 



2E P ) 2 \2gn c ) 



where E p ~ cp and c = ygn c /m is the speed of Bogoliubov sound. The sign of v v in (|8l|) is 
opposite from the one given in pO| , pT| because we have defined u p v p in fl55|) to be positive. 
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In this limit, one obtains for the scattering amplitude the following expression [19 



This approximation is valid at low temperatures, where only low-momentum excitations are 
relevant. 

3) Finally, one can consider the scattering of phonons (low-momentum excitations) with 
momentum p% -C po with particles (high-momentum excitations) with momenta Pi,P2 3> Po- 
The corresponding amplitude for this process is given by |19| 



A(l;2,3)^2 3 / 4 (p 3 M,) 1/2 . (83) 

Therefore, in the case of a sound wave scattering with particle-like excitations the scattering 
amplitude only depends on the wavevector of the sound wave. It is independent of the 
momenta of the scattering particles. 

V. LOCAL EQUILIBRIUM SOLUTION 

To describe the thermalization of quasiparticles, it is sufficient to consider the C22 colli- 
sion integral. From ([48), (|j2]) and (|62|) , one obtains 



C22U) = -^g 2 J -^^ S(E P + cji - uj 2 - uj 3 )5(p + pi - p 2 - p 3 ) 



1 2 f dpiduoi 
x [//i(l + / 2 )(1 + h) - (1 + /)(! + /i)/ 2 / 3 ] 



[Tr (a(p 2 , ^2)a(p, E p )) Tr (a(pi, Wi)a(p 3> w 3 )) 
+ Tr (2d(p 2 , u 2 )a(p 1 , Wi)d(p 3 , w 3 )d(p, £? p ))] . (84) 

The local equilibrium distribution function for quasiparticles /o(p, o;;R, T) is determined 
by the requirement that C^t/o] = 0. One can see from (^) that one doesn't have to 
specify some specific approximation for the single-particle spectral densities. We only need 
a solution for / such that the expression in fl54"| ) containing the /'s vanishes. One can verify 
that C*22[/o] — if /o(p, a;; R, T) has the following form 
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/o(p, u; R, T) = e ^_ p , (Vn _ Vs) 1 _, 9p(R , r)) _ 1 - (85) 

The vector v n — v s describes the mean drift velocity of the quasiparticle gas in the local rest 
frame of the superfluid and for small velocities we have P,|H| 

Pn (v n - v s ) = J dpp/ (p, R, T 7 ). (86) 

which defines the normal density p n . This is consistent with the usual Landau definition. 
The distribution function given in (|55D differs from the usual equilibrium quasiparticle distri- 
bution function discussed in the standard literature for phonons and rotons in liquid helium 
|14| . Since the number of quasiparticles in not conserved, the usual form for the equilibrium 



quasiparticle distribution function has no chemical potential (i.e., the chemical potential is 
zero). Here, we have introduced a chemical potential in ([85]) to allow for the possibility 
that the condensate atoms and the thermally excited quasiparticles are not in the diffusive 
equilibrium with each other. 

To understand the physics of the quasiparticle chemical potential in (|85|) better, let 
us first consider the high temperature case. At high temperatures, the particle and the 
quasiparticle excitation spectrum are equivalent and the the local equilibrium distribution 
function in the lab frame is given by || 

/o(p, R, T) = -j- , (87) 

^( '"7:"' +t/(R,r)-A(R,r)j _ 

with U(R, T) = U ext {K) + 2gn{K, T) and u - v s • p = E p (see If we transform (§7p to 

the local rest frame (where v s = 0), then p' = p — mv s is the momentum in the local rest 
frame and fl87|) becomes 

/o(p, R, T) = — = : . (88) 

Here E p > = |— + gn c (H, T) is the excitation energy in the local rest frame. Therefore, if we 
define @ 

Hdiff = ~ He - ^m(v n - v s ) 2 , (89) 
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we see the quasiparticle chemical potential fi qp introduced in ( |85D can be identified with Hdiff 
discussed at length by Zaremba, Nikuni and Griffin ||. Hence, we see that in the quasipar- 
ticle description in a local frame in which the condensate is stationary, the difference {fJ-diff) 
between the chemical potentials of the condensate and non-condensate that was introduced 
by ZNG to describe the non-diffusive equilibrium of these two components appears very nat- 
urally as the quasiparticle chemical potential. The standard case discussed in the superfluid 
helium 4 He literature fl4| corresponds to n qp = (see, however, the discussion of the second 



viscosity coefficients in superfluid 4 He fUjl ). 

To summarize, we can distinguish two distribution functions / which satisfy C22L/0] — 0: 

(1) The condensate atoms and the quasiparticle excitations are in diffusive thermal equi- 
librium, i.e. n qp = and hence Ci 2 [/] = 0. 

(2) The condensate atoms and the quasiparticle excitations are not in the diffusive thermal 
equilibrium, i.e. ji qp 7^ 0. In this case, one finds that Ci2[/] is proportional to [1 — e^ M?p ], as 
in ZNG §. 

VI. GENERALIZED GROSS-PITAEVSKII EQUATION 

As we noted in the Introduction, the dynamics of a trapped Bose-condensed gas is usually 
described in the literature by the Gross-Pitaevskii equation of motion ([!]) for the condensate 
order parameter. Linearizing this equation, one obtains the collective mode frequencies that 
have been confirmed in many experiments However, at finite temperatures, the sim- 
ple GP equation does not provide an adequate description of the thermally excited atoms. 
Moreover, even at T -C Tbec-, m recent experiments at JILA on 85 Rb f34| , the dimensionless 



gas parameter can be as large as y/n c a 3 ~ 10 _1 , i.e. the quantum depletion of the conden- 
sate is not negligible. The simplest generalization of the zero temperature GP equation of 
motion is usually done by including an additional self-consistent Hartree-Fock mean field 
2gn produced by the thermally excited atoms. The condensate atoms described by $(R, T) 
move in this mean-field, in addition to the field produced by the condensate. However it 



29 



is clear that the second-order collisions which we have included in deriving the quasiparti- 
cle kinetic equation in Section IV must also be included in a generalized GP equation for 
$(R, T). Technically, this arises from the three-field correlation function given in (p3|) . 

Zaremba, Nikuni and Griffin || have evaluated this three-field correlation function for 
the thermally excited atoms following the method of Kirkpatrick and Dorfman |5|, a method 
which is not very transparent. ZNG obtained a generalized Gross-Pitaevskii equation with a 
dissipative term associated with the Cu collisions. This new dissipative term is, as expected, 
proportional to the collision integral for scattering between atoms in the condensate and 
thermal atoms, since such collisions change the number of atoms in the condensate. The 
ZNG work was limited to finite temperatures where the thermal atoms can be described 
as free atoms moving in the dynamic Hartree-Fock mean-field produced by all other atoms 
(both those in the condensate and in the thermal cloud). Recently, we have derived the 
same generalized GP equation as ZNG using the powerful KB method method. In this 
section,we now extend this kind of calculation to deal with low temperatures. The new 
equation of motion for the order parameter will be shown to be identical to that obtained 
in Refs. [|]|| apart from the fact that C 12 is now given by the expression in ([77]). That is, 
it now involves the Bogoliubov quasiparticles and collision cross-section is renormalized by 
various Bose-coherence factors involving the w's and v 's. 

To derive an equation of motion for the condensate order parameter, we first write 
equation (|30|) for $(r, t) in the local rest frame. As before, under the gauge transformation 
(|32|), the only change is that the non-interacting propagator is now given by (|34]) . The 
equation of motion in the new local frame is (see (|27|)) 

$(1) 

= f dl [S> - S<] (r-f,t- t; (r + f)/2, (t + t)/2) $(f , t) 

J —oo 

+ f dl[S>-S<](r-T,t-i-,(r + f)/2,(t + t)/2)®*(r,t). (90) 

J — oo 

Here, we have rewritten the condensate self-energy in the center-of-mass and relative coor- 
dinates and, as usual, set in = (the Popov approximation). We recall that in the local rest 
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frame, the order parameter phase is removed and hence $(r, t) = yn c (r, £). 

We assume, as usual, that the S correlation function (PS|) is dominated by small values 

> 

of the relative space-time coordinates (r — r, t — t). Hence we can approximate in ( |9"U| ) 

> 

by S < (r — f, t — t;r,t). For the same reason, we can also approximate the macroscopic 



wavefunction $(r, i) in the integrand of fl90|) by $(r, t) = yn c (r,t). Hence, ( P0|) simplifies 
to 



f dfdt[(S> - S<) (r - r, * - i; r, t) $(r, t) + (S> - S<) (r-f,t- t; r, t) $*(r, *)] . (91) 

■/ — oo 



We can rewrite ( pT|) (labeling (r, £) — > (R, T)) as follows 



d de(R,T) 1 r _ TT ^ 

i-^ - K dT } +-[V R + ™v s (R, T)] 2 + /i - EU(R) 



s(2n(R,T) + n c (R,T))]$(R,T) 

R, T) / ^ [S> - 5< + 5> - 5<] (p, w; R, r) jT dHfe*^" 4 ^. (92) 



In the second-order Beliaev approximation, the condensate self-energy is given by ( p9f) . The 
Fourier transform of this is 



S<(p, w; R, T) = -^g 2 J d ^^ 5(u + u x - u 2 - u 3 )5(p + Pi - p 2 - p 3 ) 

■ < > 
<P(pi, Wi; R, T)^<(p 3 , cj 3 ; R, T) 



(27T) 

> 

^ < (P2,^2;R, T)Tr 
+ 2~g<(p 2 , u 2 - R, T)^>( Pl , Wl ; R, T)r (p 3 , u; 3 ; R, T) 



(93) 



In evaluating the right-hand side of (p2|) , we use the identity 

lim [ T die-^ 5 ^-^ ~ 7c5(co) + iP (-) 

S^O+J-oo \UjJ 

and only keep the delta function part, to obtain 

d (90 (R, T) 1 r _ , n _ Nl2 TT tns 

l df dT + 2^ [Vr + zmv *( R ' T )] + ^ ~ ^( R ) " 

- (5> - S< + S> - S<) (p = 0, u = 0; R, T)] $(R, T) = 0. 



(94) 



g(2n(R,T) + n c (K,T)) 
(95) 
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Using (f48|) in the condensate self-energy S given in fl93|) , the second-order terms appearing 
in (|95| ) reduce to 



(Su - S£ + S'fa - S< 2 ) (p = 0, w = 0; R, T) = i^g 2 J -^^5{u x -u) 2 - w 3 )£(Pi - P2 - Ps) 



.1 2 f dpiduji 

x + / 2 )(1 + / 3 ) - (1 + / x )/ 2 / 3 ] [(a n + a 12 )(p 2 , u 2 )Tr (a( Pl , w 1 )a(p 3 , w 3 )) 



+ 2 (a(p 2 , w 2 )a(Pi, wi)a(p 3 , a; 3 )) n + 2 (a(p 2 , w 2 )a(pi, Wi)a(Ps, w 3 )) 12 ] . (96) 



Recalling that in the local rest frame we have $(R, T) = Jn c (R, T), with no phase, we 
finally obtain a generalized Gross-Pitaevskii equation in the following form 



d9(R,T) 1 r _ _ xl2 



dy/n c (R,T) 
1 dT 

+ U ext (R) + g [2fi(R, T) + n c (R, T)] - zi?(R, T)] ^/n c (R, T). (97) 

The new dissipative term R in the GP equation is clearly related to the Cyi collision term 
in the kinetic equation fl80D , namely [||,|J 

j^r^Z* <W<**>n. ( 98) 

V ' ; J (2tt) 3 2n c (R,T) V ; 

This term describes the damping of condensate amplitude fluctuations due to collisions with 
the thermal excitations. The appearance of the dissipative term in ( |97|) is expected since 
the C\2 collisions change the number of atoms in the condensate and hence can modify the 
magnitude of the condensate macroscopic wavefunction. We note that since we ignore the 
real part of the second-order self-energies, the condensate chemical potential in ( fj3|) is not 
modified. If we transform back into the lab frame (where we have $ = ^/n^e 10 ), ( p7|) reduces 
to the time-dependent generalized Gross-Pitaevskii equation for $(R, T) discussed by ZNG 
||. However, Ci 2 now involves the Bogoliubov quasiparticle spectrum in place of the HF 
particle-like spectrum used in Ref. ||, and in addition, the collision integral matrix elements 
involve the characteristic Bose coherence factors u and v. 
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VII. KOHN MODE 



In this section, we show that the non-condensate and condensate both exhibit the rigid 
in-phase oscillations, the Kohn mode. This mode is discussed in detail in Section VI of [Q, 
and the analysis there is easily generalized to the more general equations we are discussing. 
The center-of-mass oscillation of the non-condensate and condensate density profiles corre- 
sponding to the Kohn mode is given by 

n c (R,T) =n c0 (R-77(T)) 

n(R,T) =n (R-T7(T)). (99) 

Here, the center-of-mass displacement rf{T) (with v s = if) satisfies the harmonic oscillator 
equation of motion 

^ Va 2 / 1m \ 

m ~&r I = ~ UaT]a ' ( 10 °) 

where u a is the trap frequency in the a th direction. The quasiparticle distribution function 
/(p, R, T) corresponds to the equilibrium density profile oscillating around its center of mass 
with the trap frequency, i.e., 

/(p,R,r) = / (p,R-T7(T)). (101) 

To prove ([TOll) , we note that with ©, the expression for the Bogoliubov excitation 
energy in ([57]) reduces to 



E P (R, T) = J el + 2sn c0 (R - r))e p = E p0 (R - rj). (102) 



Therefore, the kinetic equation for the quasiparticle distribution function in ([80]) is 



^ + V p (E p0 (R - n{Tj) + r) • P) • V R - V R E p o(R - r/(T)) • V p ] / (p, R - r/(T)) 
= C 12 [/ (p,R-r;(T))]. (103) 

If we expand /o(p, R — around rj = 0, 
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/o(p, R - ry(r)) = / (p, R) - r} ■ Vr/ (p, R), (104) 
and neglect the quadratic terms in 77, (|103|) simplifies to 

V p E p0 (R - rj(T)) ■ Va/o(p, R - rj(T)) - V R E p0 (R - ry(T)) • V p / (p, R - ry(T)) 

= C la [/o(p,R-f;(T))]. (105) 

The left-hand side of ( |105|) is seen to be the kinetic equation for the equilibrium distribution 
function. To prove that the Kohn mode is a solution, one only has to show Ci2[/o(p, R — 
T](T)] = 0. Assuming the equilibrium quasiparticle distribution function f is given by (|85"D 
with \i qv = 0, and using the identity for the Bose distribution function 

l + f(x) = e x f(x), (106) 

one obtains following expression for the C 12 [/o] 

C l2 [f (p, R - „(T))] = 2^V(R - rj(T)) J dp ff P3 | A (2, 3; 1) | 2 
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x5(pi - p 2 - P3 )5 [E 10 (R - ry(T)) - ^ 20 (R - ry(T)) - ^ 3 o(R - r)(T))} 

X [5(p - Pi) - Sip - p 2 ) - S(p - p 3 )] (1 + Ao)/20/30 

x J _ e -/3[-Blo--B20--B 3 o-(pi-p2-p3)(Vn-V s )] 



(107) 



Using the delta functions in ( J107 ) corresponding to the conservation of energy and momen- 



tum, it immediately follows that 

C 12 [/ (p,R-77(T))] = 0. (108) 

This proves that the non-condensate exhibits a rigid simple harmonic displacement with the 
trap frequency. 

Since we have proven that the collision integral vanishes for the Kohn mode type of 
oscillations, the dissipative term R in the generalized GP equation in (p7| ) will vanish as 
well. Therefore, it follows that the equilibrium condensate profiles oscillate with the trap 
frequency (for the analogous calculation based on the HF and HFB single-particle spectrum, 
see Refs. 
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VIII. CONCLUSIONS 



Using the Kadanoff-Baym non-equilibrium Green's function formalism [11,12] we have 
derived a kinetic equation for the quasiparticle distribution function and a generalized Gross- 
Pitaevskii equation valid at all temperatures. Our new kinetic equation involves the Bo- 
goliubov quasiparticle spectrum and Bose coherence factors involving the u and v functions. 
Our equations reduce to those obtained at high temperatures by ZNG 0. As we have em- 
phasized in Section II, the approximation that we have used in this paper is gapless and it 
gives the correct low-momentum (long- wavelength) limit. In contrast, the so-called "con- 
serving approximations" are based on a functional from which both self-energy £ and the 



source r/ functions can be derived by functional differentiation (see p. 338 ff of Ref . [ 25[] , and 
also Section III in Ref. |4]]) . The resulting single-particle Green's function can be used to 
generate a density response function whose spectrum is guaranteed to satisfy conservation 



laws [|i8| , p5] , |35| , |36| , even though the generating Green's functions have an energy gap in the 
long- wavelength limit. 

The simple quasiparticle approximation that we have used in this paper has allowed 
us to derive the kinetic equation in a "Boltzmann-like form". The kinetic equation in the 



Bogoliubov-Popov approximation given in fl80| ) is only valid in this quasiparticle approx- 
imation. In deriving these results, we have neglected the real part of the second-order 
self-energies that give rise to many-body corrections. Our discussion could be generalized to 
include these real parts, but this improved theory would be very complex. A first step would 
be to include such renormalization effects within a simple quasiparticle approximation to 
the spectral densities a af3 . 

Our kinetic equation for the quasiparticle distribution and the generalized Gross- 
Pitaevskii equation are coupled and have to be solved self-consistently. They should provide 
a sound basis for the future systematic study of the non-equilibrium response of a trapped 
Bose gas at low temperatures. We remark that we could use our results to derive the 
Landau-Khalatnikov two-fluid hydrodynamic equations, in the collision-dominated region. 
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Indeed, the approach developed in Ref. [14] is based on a quasiparticle kinetic equation 



which is precisely of the kind we have derived in (|50|). Such a calculation would extend a re- 
cent derivation |TI] of the Landau-Khalatnikov two-fluid equations in the high temperature 
region where the simple HF single-particle spectrum is appropriate. 
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APPENDIX A: 

For illustration, we give the detailed proof that the collision integral / in (^) conserves 
momentum, namely 

J dpp/[/(p,R,T)] = 0. (Al) 



The proof is essentially the same as one uses in classical gases (See Ch. 5 of [0). Using the 



expression for the second-order Beliaev energy in (|S2|) in conjunction with (fig ), one obtains 



f , rr _, 1 2 f dpdu dpidui 

J dppl[f{p, R, T)\ = --g J 5{u + u x - u 2 - u 3 )8(p + p x - p 2 - p 3 ) 

P [[//i(l + / 2 ) - (1 + /)(1 + /O/a/a] [Tr (a(p 2 , u 2 )a(p, u)) Tr (a( Pl , Wl )a(p 3 , w 3 )) + 

2Tr (a(p 2 , u 2 )a(pi,u 1 )a(p 3 , uj 3 )d(p, u))) + (A2) 

[f(l + / 2 )(1 + / 3 ) - (1 + f)hh] [Tr (a(p 2 , cu 2 )a(p, u)) Tr (h(pt, Wi)a(p 3 , w 3 )) + 

2Tr (a(p 2 ,a;2)^(pi,a;i)a(p 3 ,a; 3 )a(p,a;)a(p,a;)J + (A3) 

[//i(l + / 2 ) - (1 + /)(1 + fi)f 2 ) [Tr (a(p 2 , a; 2 )a(p, u)) Tr (a( Pl , ^)Mp 3 , w 3 )) + 

2Tr (a(p 2 ,uj 2 )a(p 1 ,uj 1 )h(p 3 ,u 3 )a(p 1 u)a(p,uj)j + (A4) 
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[//i(l + /s) - (1 + /)(1 + /i)/ 3 ] [Tr (A(p 2 , W2)a(p, w)) Tr (a( Pl , c*)a(p 3 , w 3 )) + 
2Tr (h( 



p 2 , cj 2 Ja(Pi, WiJa(P3, w 3 )a(p, v)a(p, u > 



Consider the ( |A2|) term first. The change of dummy variables 



(A5) 



(P3,^3j 



(A6) 



doesn't change the delta functions. Because of the cyclic invariance of trace, the first term 



in (|A2|) is also unchanged - the only change is that p in front of the integral becomes — p 2 . 
After the change of variables given by (|A6|), the second term in (|A2|) becomes 



Tr (a(p, w)a(p 3 , ct» 3 )a(pi, Wi)a(p 2 , u; 2 )) = Tr (a(p 2 , c^ 2 )a(p, w)a(p 3 , u; 3 )a(pi, a>i)) (A7) 



Transforming (p, cj) — ► (pi, a>i) in (|A7|), doesn't change either the 5-functions or the product 
of /'s. Therefore the trace in second term in ( |A2| ) doesn't change. Thus (|A2|) after the 
transformations is unchanged, but p is replaced with — p 2 . Using the transformation, 



v Pi,o;ij (p 2 ,o; 2 J 



vP3,^3j 



(A8) 



we would again obtain the same expression as (|A2), but with p replaced with p x . Finally, 
if we make the transformation 



(p, a;) ^ (p 3 , w 3 ) (P2, w 2 ) ^ (pi, u)\) 



(A9) 



we obtain the same expression as in (|A2|) , but with p replaced with — p 3 . We conclude that 
we can write ( |A2|) as one fourth of the sum of four equivalent terms. The integrand of this 
new expression is thus seen to be proportional to 



(P + Pi - P2 - Pa) 8 (p + Pi - p 2 - p 3 ) 

which clearly vanishes. 

A similar discussion can be given of the other terms in 
of variables 



(A10) 



). Making the change 
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(p,w) ^ (pijWi) (p 2 ,^> 2 ) ^ (p 3 ,^ 3 ) (All) 



in ( |A4] ) first, and using the cyclic invariance of the trace, one can show that ( |AJ ) is the same 
as (|A5|) , but with p x instead of p as an overall multiplying factor. Using 

(p,w) ^ (P3,W 3 ) (P2,W 2 ) ^ (P2,^2)- (A12) 

in (|A3|) , and then relabeling p — > pi in the second term in (|A3|) , we can reduce ( \K3j ) to the 
same expression as in ( |A5| ) but with p replaced with — P3. Therefore the sum of the terms 
), (|A4 ) and (|K5|) has an integrand involving 



(P + Pi - Ps) 5 (P + Pi - P2 - P3) <5(p 2 )5(cj 2 ), (A13) 

which clearly vanishes. This completes the proof of (|A1|) and hence of momentum conser- 
vation by collisions. 



38 



REFERENCES 

[1] F. Dalfovo, S. Giorgini, L.P. Pitaevskii, and S. Stringari, Rev. Mod. Phys. 71, 463 
(1999). 

[2] See, for example, P. Sokol, in Bose-Einsten Condensation, ed. by A. Griffin, D.W. Snoke 
and S. Stringari (Cambridge University Press, N.Y. 1995). 

[3] E. Zaremba, T. Nikuni, and A. Griffin, Journ. Low Temp. Phys. 116, 277 (1999); for a 
brief account of this work, see also T. Nikuni, E. Zaremba and A. Griffin, Phys. Rev. 
Lett. 83, 10 (1999). 

[4] M. Imamovic-Tomasovic and A. Griffin, in Progress in Nonequilibrium Green' Functions, 
ed. by M. Bonitz (World Scientific, Singapore, 2000), p. 404; e-print |cond-mat/9911402 . 

[5] H.T.C. Stoof, Journ. Low Temp. Phys. 114, 11 (1999). 

[6] R. Walser, J. Williams, J. Cooper and M. Holland, Phys. Rev. A 59, 3878 (1999). 

[7] C.W. Gardiner and P. Zoller, Phys. Rev. A 61, 033601 (2000). References to earlier 
work by these authors are given here. 

[8] N.P Proukakis, K. Burnett and H. T. C. Stoof, Phys. Rev. A 57, 1230 (1998). 

[9] T.R. Kirkpatrick and J.R. Dorfman, Journ. Low Temp. Phys. 58, 308 (1985); 58, 399 
(1985). 

[10] L.P. Kadanoff and G. Baym, Quantum Statistical Mechanics (W.A. Benjamin, New 
York, 1962). 

[11] Hua Shi and A. Griffin, Phys. Reports 304, 1 (1998). 

[12] J.W. Kane and L. P. Kadanoff, Journ. Math. Phys. 6, 1902 (1965). 

[13] J.W. Kane, Ph.D. Thesis, University of Illinois, Urbana, 1966. 

[14] I.M. Khalatnikov, An Introduction to the Theory of Superfluidity (W.A. Benjamin, N.Y., 

39 



1965). 

[15] D. Zubarev, V. Morozov and G. Ropke, Statistical Mechanics of Nonequilibrium Pro- 
cesses, (Akademie Verlag, Berlin, 1997), Vol. II, Ch. 6. 



[16 

[17; 

[18 
[19 
[20 
[21 
[22 
[23 

[24 
[25 
[26 
[27 

[28 

[29 
[30 
[31 



S. Giorgini, Phys. Rev. A 61, 063615 (2000). 

P. O. Fedichev and G. V. Shlyapnikov, Phys. Rev. A 58, 3146 (1998). 

A. Griffin, Phys. Rev. B 53, 9341 (1996). 

U. Eckern, Journ. Low Temp. Phys. 54, 333 (1984). 



M.J. Bijlsma, E. Zaremba and H.T.C. Stoof, Phys. Rev. A, in press; |cond-mat/0001323| 2 



T. Nikuni and A. Griffin, submitted to Phy. Rev. A, |cond-mat / 0009333 



M. Imamovic-Tomasovic and A. Griffin, Phys. Rev. A 60, 494 (1999). 

A.L. Fetter and J.D. Walecka, Quantum Theory of Many Particle Systems (McGraw- 
Hill, New York, 1971). 

S.T Beliaev, Soviet Phys. JETP 7, 289 (1958). 

PC. Hohenberg and PC. Martin, Ann. Phys. (N.Y.) 34, 291 (1965). 
T.H. Cheung and A. Griffin, Phys. Rev. A 4, 237 (1971). 

N.N. Bogoliubov, Lectures on Quantum Statistics (Gordon and Breach, N.Y. , 1970), 
Vol. 2. 

H. Haug and A. -P. Jauho, Quantum Kinetics in Transport and Optics of Semiconductors 
(Springer, Berlin, 1996), Part II. 

G. D. Mahan, Many-Particle Physics (Plenum Press, 1990), 2nd ed., Sect. 7.4. 

S. Giorgini, Phys. Rev. A 57, 2949 (1998). 

L.P Pitaevskii and S. Stringari, Phys. Lett. A 235, 398 (1997). 

40 



[32] A. A. Abrikosov, L.P. Gorkov and I.E. Dzyaloshinski, Methods of Quantum Field Theory 
in Statistical Physics (Dover, N.Y., 1963), p. 11 ff. 

[33] E. M. Lifshitz and L.P. Pitaevskii, Statistical Physics, Part 2 (Butterworth and Heine- 
mann, Oxford, 1995), Ch. 3. 

[34] S.L. Cornish, N. R. Claussen, J.L. Roberts, E.A. Cornell and C.E. Wieman, Phys. Rev. 
Lett. 85, 1795 (2000). 

[35] G. Baym and L.P. Kadanoff, Phys. Rev. 124, 287 (1961). 

[36] G. Baym, Phys. Rev. 127, 1391 (1962). 

[37] K. Huang, Statistical Mechanics, 2nd ed. (J. Wiley, New York, 1987). 



41 



